Heat transport by laminar boundary layer flow with polymers 
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Motivated by recent experimental observations, we consider a steady-state Prandtl-Blasius bound- 
ary layer flow with polymers above a slightly heated horizontal plate and study how the heat trans- 
port might be affected by the polymers. We discuss how a set of equations can be derived for the 
problem and how these equations can be solved numerically by an iterative scheme. By carrying out 
such a scheme, we find that the effect of the polymers is equivalent to producing a space-dependent 
effective viscosity that first increases from the zero-shear value at the plate then decreases rapidly 
back to the zero-shear value far from the plate. We further show that such an effective viscosity 
leads to an enhancement in the drag, which in turn leads to a reduction in heat transport. 

PACS numbers: 



I. INTRODUCTION 



It has been known for more than 60 years that adding 
polymers into turbulent wall-bounded flows can reduce 
the friction drag significantly (see, for example, [H, 0] 
and references therein). This effect of polymer additives 
on mass transport has been studied extensively during 
the past 60 or so years. On the other hand, the effect of 
polymers additives on heat transport is much less stud- 
ied. Recently, an experimental study reported Q that 
for turbulent Raylcigh-Bcnard (RB) convection of wa- 
ter, confined within a cylindrical cell heated below and 
cooled on top, adding polymers to the flow reduces the 
heat transport. In turbulent RB convection within a con- 
tainer of given shape, the flow state is characterized by 
two dimensionless parameters: the Rayleigh number (Ra) 
and the Prandtl number (Pr), where Ra measures the size 
of the thermal forcing and Pr is the ratio of the kinematic 
viscosity to the thermal diffusivity of the fluid. Moreover, 
there is an exact balance Q between the heat transport 
and the energy and thermal dissipation rates. The en- 
ergy and thermal dissipation rates can be decomposed as 
sums of contributions from the bulk of the flow and from 
the boundary layers near the top and bottom plates Q. 
The experimental study reported Q was conducted at a 
Ra of the order of 10 10 . At such a moderate Ra, the en- 
ergy and thermal dissipation rates are mostly due to the 
contribution from the boundary layers @ . This suggests 
that the observed reduction in heat transport is likely to 
be an effect of the polymers interacting with the bound- 
ary layer flow in turbulent RB convection. Experimental 
measurements @ indicated that the average velocity and 
temperature boundary layer profiles in turbulent RB con- 
vection at moderate Ra could be described by those pro- 
files in the steady-state Prandtl-Blasius boundary layer 
flow @, II] above a slightly heated flat plate. 

Motivated by these experimental observations, we 
study a steady-state Prandtl-Blasius boundary layer flow 



with polymers near a slightly heated plate and focussed 
particularly on the possible effect of the polymer addi- 
tives on the heat transport. Physically, we can think 
of this boundary layer flow as the flow near the bottom 
plate of the convection cell in turbulent RB convection. 

This paper is organized as follows. In Section [TTl we 
formulate the problem and set up the equations of mo- 
tion for the system. We discuss how the set of equations 
can be solved numerically by an iterative scheme in Sec- 
tion Mil After solving the problem using such a scheme, 
our results show that the effect of the polymers is equiv- 
alent to producing a space-dependent effective viscosity 
that first increases from the zero-shear value at the plate 
then decreases rapidly back to the zero-shear value far 
from the plate. We further show that such an effective 
viscosity would lead to an enhancement in the friction 
drag and a reduction in heat transport. We shall present 
and discuss our results in Section IIVI Finally, we shall 
summarize and conclude in Section FVl 



II. THE PROBLEM 

For the Prandtl-Blasius boundary layer flow above a 
large flat plate, the velocity equation is: 



v x d x v x + v v d v v x 



vd yy v x 



(1) 



Here x denotes the direction along the plate, y denotes 
the direction away from the plate, and v is the kinematic 
viscosity of the fluid. Both v x and v y vanish at the plate 
and v x merges with the uniform mainstream velocity U 
far away from the plate. The crucial point about Eq. ([1]) 
is that the viscous term is balanced against the nonlinear 
advection term and that the flow changes more rapidly 
away from the plate than along the plate such that d y 3> 
d x . The latter condition is satisfied for flows with large 
Reynolds number. 
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Introducing the variable 

£ = \/ — v 

V vx 

and the stream function 

^(x,y) = \/vxU<t>(£) 

such that 



(2) 



(3) 



(4) 



(5) 



showing that the velocity profile is self-similar in the di- 
mensionless variable £ at different position x along the 
plate. Here denotes d^<p. The boundary conditions 
are 



d x ^> 



we obtain the famous Blasius equation Q 

20«£ + <M>& = o 



0(0) = e (O) = ; f (oo) = l 



(6) 



as v x — > U when it is far away from the plate. The plate 
is slightly heated at a temperature T\ above the ambient 
temperature To far away from the plate. Writing the 
temperature field as 



T(x,y)=T + (T 1 -T )9(£) 



(7) 



[Relating to turbulent RB convection, To can be taken 
as the temperature at the center of the cell and thus the 
total temperature difference across the cell is 2(Ti — To)]. 
Then 8 satisfies the equation 



2% + #> e Pr = 



where Pr = v/k and n is the thermal diffusivity of the 
fluid, and the boundary conditions are: 



0(0) = 1 



(oo) = 



(9) 



We want to investigate the effect of polymers on heat 
transport in this laminar Prandtl-Blasius flow. The poly- 
mers produces an additional stress in the momentum 
equation of the fluid. This polymer stress 7y depends 
on the amount of stretching of the polymers and is thus 
a function of the dimensionless conformation tensor Rij 
of the polymers. Let the vector d denote the polymer 
end-to-end distance and po be the polymer radius in the 
unstretched regime, then represents the average over 
many (TV) polymers in a small region around the point 
(x, y) of the product didj/p^, i.e. Rij = N~ 1 T,didj / p\. 
In the simplest Oldroyd-B model of polymers [|[ , 



7- (Ri. 



Sij) 



(10) 



where v p is the polymer contribution to the viscosity of 
the solution at zero shear and r is the relaxation time 
of the polymers. Thus in the presence of polymers, the 



equation of motion for the velocity field is modified by 
an additional stress that depends on R^. By employing 
the same ideas leading to the Blasius equation [Eq. |(T}] 
we have 



v x d x v x + v v d v v x 



vd yy v x 



On 



P p 



(11) 



It has been shown that the relaxation time of the poly- 
mers can be significantly increased by the stretching of 
the polymers [1(|. To model this effect, we let 



r = r 



1 + aR 
1 + a 



(12) 



where To is the bare Zimm relaxation time, R = (R xx + 
Ryy) 1 / 2 , and a > is a parameter. When there is no 
stretching, R = 1 and r reduces back to To. When there 
is stretching, R > 1 and r > tq. 

In the presence of polymers, the transformation using 
£ docs not lead to a similarity solution in general in that 
explicit appearance of x remains in the equation for 0. 
This is known in the literature. Similarity solution has 
been obtained in some special cases with certain velocity 
or temperature boundary conditions (Til IT2I ] that might 
not have direct physical relevance. Here to circumvent 
this difficulty, we recall that the Prandtl-Blasius flow is 
meant to be applicable when x is large (such that d y ^> 
d x ). Thus we make the following approximations: 



L 



m - 0) 



(13) 



—LA 
° x ~ 2Td£ 



(14) 



(8) by putting x = L, the length of the (long) plate, and 



replacing v in £ [Eq. (|2J)] by vq 



where v$ is 



the total viscosity of the polymer solution at zero shear. 
That is, in the presence of polymers, we have 



(15) 



£ = \ y with polymers 

vox 



With these approximations, the scaling transformation of 
£ leads to a similarity solution. The resulting modified 
Blasius equation is: 



^ = (l- 7) ^ e + ^=| 



{l + a)R x 



aR 



(16) 



where the Weissenberg number (Wi) and the Reynolds 
number (Re) are defined as 



Wi 



t U 
L 



Re 



UL 



(17) 



and 7 = fp/fQ is a function of the polymer concentration. 
As usual in the Prandtl-Blasius approximation, all terms 
of the order of 1/Re are nelgected in Eq. p^|) . 

We would like to study how the heat transport is af- 
fected by the polymers. In turbulent RB convection, it 
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is common to measure the heat flux Q in terms of the 
dimensionless Nusselt number (Nu) , which is the ratio of 
Q to that when there is only conduction, defined by 



Nu 



Q 



dT 
dy 



y=0 



2fc(Tj - T )/H 2(Tj - T )/H 



(18) 



where k is the thermal conductivity of the fluid, H is the 
height of the convection cell, and (. . ) a is the average 
over the cross section of the cell. For the Prandtl-Blasius 
flow, taking H — L and dropping the numerical factor, 
Nu can be estimated as 



Nu 



■(0)] 



(19) 



To proceed, we must supplement Eq. (1161) with a spe- 
cific information on R xy . In a fluid flow of velocity v, 
the components of the dimensionless polymer end-to-end 
distance U = di/po, i = x,y obey the differential equa- 
tions: 



dh 



1 



= ——(h — loi) + IjdjVi + thermal noise (20) 
at It 

where Iq x = cos a and lo y — sin a and a is a random angle 
uniformly distributed in [0, 27r]. Neglecting the thermal 
noise, we can rewrite Eq. pOj) for the two-dimensional 
Prandtl-Blasius flow as 



L dl x 
U~dT 
L dl y 
U~dt 



(1 + a)(l x - l 0x ) 

2Wi(l + aR) 2 
(1 + a)(l y - l 0y ) 



2Wi(l + aR) 



4VRe 



Re<^(21) 
%,(22) 



In order to obtain U as a function of £, we assume that 
each polymer follows the streamline of the flow such that 
d£/dt = v x d x ^ + Vydy^ = —(U/2L)<f>, again using the 
approximations in Eqs. (|lUp and f|14|) . Thus (L/U)d/dt = 
-(cj)/2)d/d4 and Eqs. (J2TD and ^ become 



dl x 

(My 



(1 + a)(l x - l 0x ) 
Wi(l + aR) 

(I + a){ly - l y) e<Ptt 

Wi(l + aR) 2\/Re" 



^l x - 2vTW> 5? Z y (23) 

L - Mttly (24) 



The quantity R xy is the average of l x l v over all the 
polymers contained in a small volume centered near the 
point (x, y), over the angle a, and over the thermal noise. 
Such an average would depend on the precise distribution 
of the polymers, which is not obvious to obtain. Thus in- 
stead of performing such an average, we calculate l x and 
l y for one polymer for a fixed angle a = 7r/4 such that 
lox = k)y = l-o = l/v2) and estimate R xy as some func- 
tion of the calculated l x and l y . In an earlier short com- 
munication [i~3| , we have shown that neglecting the ther- 
mal noise and in the limit of small Wi and a = 0, the av- 
eraging over the angle a gives R xy = [l + /i(£)]</>££Wi\/Re 



with h(£) — >• as £ — >• oo. We expect this behavior of R xy 
going like Wiv / Re</>££ far away from the plate is gener- 
ically true. Moreover, far away from the plate, R — >■ 1 
and l y — > Iq. Thus we estimate R xy as 

-^)ct>tt + (l x -l )(ly-lo) (25) 



R xy = WiVRc 



1 



1 



which has the desired asymptotic behavior far away from 
the plate. Defining 



MO-h][ly{O-k]0- + a) 



' 9(e) WiVRc[l + aR(0}^{0 
with R(0 = [l 2 x (0+l 2 y(0] 1/2 >then 

(1 + aR) 



R 



WiVfe W fe(1+fl) 

Substituting Eq. ([27]) into Eq. ([16]). we have 



27^7 







(26) 



(27) 



(28) 



Comparing Eq. (|28p with Eq. ([3]), it is obvious that 
the effect of the polymers is contained in the term with 
g(£). Furthermore, this effect of the polymers is equiv- 
alent to producing a space-dependent viscosity v s{Q- 
Spechcally, if we define v c h by 

-^-R xy (0 = Ves(0dyV* (29) 



and take 



I/eff(0 = "p[l + $(0] 



(30) 



p 1 

then we obtain exactly Eq. ([28]) from Eq. (fTTj) . The total 
effective viscosity of the polymer solution is thus 

^tot(0 = V + Ves(0 = i>0 + Vp9(Q (31) 

with the nontrivial effect of the polymers contained in 
the term Vpgiji), which is a function of the polymer con- 
centration. 

Equations $5$, (JUJ), and dHJ) should be solved 
consistently. With the obtained </>(£), Eq. ([8]) is solved 
to obtain and from which Nu is obtained by us- 
ing Eq. (fT9|). We would like to compare this Nu in the 
presence of polymers to the reference value Nuo for a 
Newtonian fluid with the same kinematic viscosity vq at 
the plate. To get the reference value Nuo, we start with 
the Prandtl-Blasius velocity profile for a Newtonian fluid 
of kinematic viscosity v$, denoted as (ft^HO- This is 
just the solution to Eq. with £ given by Eq. (TT5|) or 
Eq. (JUD with 7 = 0. With this we solve Eq. © 

and obtain the resulting temperature profile at 
the same Pr = uq/k as the polymer solution. Then 




Nu = 

We are interested in the ratio 

Nu _<%(0) 
Nuo ~ 



f\0) 



(32) 



(33) 
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lution for l x (0) is: 
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1 l 1 l l^-j 


_L 1 1 



l + Jj^>2Vfefe(0)i (35) 



As a result, for a given value of Re, the range of allowed 
values of Wi is smaller for larger values of a. Moreover, 
the allowed range of Wi is larger for a smaller value of 
Re. We have used two values of Re: 100 and 4900. For 
Re=4900, we take a = 0.01 and find that for 7 = 0.2, the 
maximum allowed value of Wi is about 2.8. For Re=100, 
we take a = 0.07 and study different values of Wi up to 
3.0 for the same value of 7. We then fix Wi=2.5, and 
study the effect of polymer concentration by varying the 
value of 7 for Re=4900 and a = 0.01. We study also 
a few other values of the parameters to investigate the 
dependence of the results on the parameters. 

IV. RESULTS AND DISCUSSIONS 



FIG. 1: Fast convergence of the numerical iterative procedure. 
Result of g(£) obtained after the first (dashed), second (dot- 
dashed), fifth (dotted) and 10th iteration (solid). 



III. CALCULATIONS 



We solve Eqs. <H2]), (HJ), and (J2SJ) consistently 
by iteration at fixed values of Re and a. We start with 
4> = in Eqs. ([23]) and ([24]) and solve for l x (£) and 



ly(£). Using the boundary condition that l x and l y 



In 



as £ — > 00, we put l x = l y = Iq at some large value 
of £, denoted as (we use £oo = 30) and integrate 
backwardly until £ equals to some finite £0 close to zero. 
We cannot integrate forwardly in £ because </> vanishes 
at £ = [sec Eq. ©]. With the calculated l x and l y , 
we obtain g(£) for £ between £ and £oo using Eq. ([26]) . 
From Eq. ([2"I]). we see that l y (0) = l and thus g(0) = 0. 
Between £ = and £0 we extrapolate g(£) using a poly- 
nomial fit. We input this g into Eq. (|28[) to solve for an 
updated 4>- Then we use this updated <p hi Eqs. (|23[) and 
([2~4"]) to obtain an updated g. Wc repeat the procedure 
until convergence in both g and <f> is achieved. 

In Fig. [T] we show g(£) obtained using this iterative 
procedure at Re= 4900, a = 0.01, Wi = 2.8 and 7 = 0.2. 
Convergence is fast and achieved after only a few itera- 
tions. Using the converged <f>, we obtain Nu for Pr=4.4 
as discussed in Section [III an d study the ratio Nu/Nuo 
for different values of the parameters. 

There is a constraint on the possible values of the pa- 
rameters imposed by the realizability of ^(0). From 
Eq. ([2"3"]) . we see that ^(0) satisfies the following equa- 
tion: 

1 + [f g (0) - h] = A(0) 2 + 11 + - (34) 



2WiVRe<^(0); 
Thus the condition for Eq. ([M]) to have a finite real so- 



We find that g(£) increases from zero at the plate 
(£ = 0) up to a certain maximum then decreases rather 
rapidly back to zero far away from the plate (£ — > 00). 
Whenever is larger than zero, the viscosity of the 
polymer solution is enhanced compared to that of the 
solvent. We expect that such an increase in the viscosity 
gives rise to an enhancement of the friction drag, which 
would in turn result in a reduction of the horizontal ve- 
locity. Indeed, we find that the horizonal velocity, v%. 
for the flow with polymers is reduced compared to v x , 
the horizontal velocity for the flow without polymers (see 
Fig. 12]). 

We study the change in drag by measuring directly the 
drag coefficient, defined by: 



C = 



y=o 



1 



l/2t/ 2 2^ 
We are interested in the ratio 



c_ 



4f (0) 



(36) 



(37) 



An enhancement in the friction drag implies a reduction 
of heat transport. To see this, we note that upon double 
integration of Eq. flSJ) by £, we obtain: 



*«(0) = 



1 



Jo dtexp[-^- J Q (j)(s)d 



(38) 



which tells us that Nu is a functional of <j>. Define $(£) = 
J Q 4>(s)ds, we can calculate <5Nu, the variation in Nu due 



to a variation in $: 
Pr 



(5Nu 



2VRc 



poo 

Nu 2 / cxp(-Pr$(s)/2)<5$(s)ds (39) 
Jo 



In the above expression, (5$ represents the variation of $ 
due to the effect of the polymers. Since v x = U<j>£, the 



5 



0.6 



<y< 0.4 



-o 0.2 
0.0 



1 1 


1 






7\ 




\ 




\ 

\ 




■f- \ 

\ 




I \ 








1 \ 




1 \ 




1 \ 
1 \ 




L 








1.1,1. 



4 

5 




8 O.Cg 



FIG. 2: Effective viscosity due to the polymers as represented 
by 73(C) (long dashed curve) and the difference in the horizon- 
tal velocity of the polymer solution from that of the pure sol- 
vent, {v£-v°)/U (solid curve), as a function of £ for Re=4900, 
a = 0.01, Wi=0.28, and 7 = 0.2. 



FIG. 4: Dependence of the effective viscosity as measured 
by 73(C) on Wi at fixed polymer concentration for Re=4900, 
a=0.01 and 7 = 0.2. From bottom to top, Wi increases from 
0.5, 1.0, 1.5, 2.0. 2.5 to 2.8. 
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FIG. 3: C/Co (squares) and Nu/Nuo (circles) as a function 
of Wi for Re=4900, a=0.01 and 7 = 0.2. 



mass throughput in the x direction across a distance £ is 
given by U </>(£), and thus a reduction in mass throughput 
implies an <5$ < 0. Therefore, Eq. (|39]) shows that drag 
enhancement, i.e., reduction in the mass throughput im- 
plies a reduction of Nu. It is indeed found that C/Co > 1 
while Nu/Nuo < 1 (see Fig.©. 

As shown in Fig. [3l the amount of drag enhancement 
and heat reduction increases with Wi. This is understood 
as the result of an increase in the effective viscosity with 



Wi. In Fig. 21 we show the dependence of g(£) on Wi. It 
can be seen that the effective viscosity increases with Wi 
while the region in which the polymers are active (i.e., 
> 0) is approximately independent of Wi. 

In Figs.[5]and[nj we show the dependence of the amount 
of drag enhancement (DE), C/Co — E an d heat reduction 
(HR), 1 — Nu/Nuo, as a function of Wi at a fixed polymer 
concentration for the two different sets of values of Re and 
a studied. It can be seen that both effects are relatively 
modest. The effect is generally larger when Re is larger. 
Both %DE and %HR increases as Wi increases and for 
the range of Wi studied, they increase quadratically with 
Wi. Interestingly, for %DE the rate of increase decreases 
with Wi and this leads to the possibility of the saturation 
of the effect at large Wi. On the other hand, for %HR, 
the rate of increase increases with Wi thus the effect on 
heat reduction is larger than that but increases for %HR. 

In Fig.0 we fix Rc=4900, a= 0.01, Wi=2.5 and show 
the extent of drag enhancement and heat reduction as a 
function of polymer concentration specified by 7. We see 
that the effect increases with 7 as expected. Moreover, 
we see a possible saturation of the effect in the limit of 
the large polymer concentration. In Fig. [SJ we show the 
dependence of 7g(£) on 7. Again we see that 75 increases 
with 7. It is this increase of the effective viscosity of 
the polymers with polymer concentration that leads to 
the increase in the extent of drag enhancement and heat 
reduction when polymer concentration increases. 

We summarize the results for %DE and %HR obtained 
for different values of the parameters in Table 1 . We see 
that the effect generally increases with a, Re, Wi and 7 
while the other parameters are kept fixed. We note that 
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0.5 1 1.5 2 2.5 3 3.5 



FIG. 5: Percentage of drag enhancement (% DE) (squares) 
and heat reduction (%HR) (circles) as a function of Wi at 
fixed 7 = 0.2 for Re=4900 and o=0.01. 




0.5 1 1.5 2 2.5 3 3.5 

Wi 



FIG. 6: %DE (squares) and %HR (circles) as a function of 
Wi at fixed 7 = 0.2 for Re=100 and a = 0.07. 

drag enhancement and heat reduction are found even for 
a = although the effect is relatively modest with only 
a few percentage of HR for 7 = 0.5 and Wi=2.5. The 
effect increases with a and interestingly the increase in 
%HR is larger than that the increase in %DE when a is 
increased. 




FIG. 7: Dependence of the effect on polymer concentration: 
%DE (squares) and %HR (circles) as a function of 7 at fixed 
Wi=2.5. 




FIG. 8: Dependence of 7<?(£) on 7 at fixed Wi=2.5. From 
bottom to top, 7 increases from 0.1 to 0.5 in steps of 0.1. 



V. SUMMARY AND CONCLUSIONS 

In this paper, we have studied the problem of heat 
transport in steady-state Prandtl-Blasius flow with poly- 
mers near a slightly heated plate. We have shown how 
a set of equations can be written for the problem and 
how this set of equations can be solved numerically by 
an iterative procedure. Our results demonstrate that the 
physical effect of the polymers is equivalent to produc- 
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ing a space-dependent effective viscosity, which increases 
near the plate. Because of this increase in viscosity, 
drag is enhanced. We have shown that such a drag en- 
hancement then leads to a reduction in heat transport 
in the Prandtl-Blasius flow with polymers. As discussed 
in Sec. Q] in turbulent RB convection of Ra about 10 10 , 
heat transport is dominated by contributions from the 
boundary layers and the mean velocity and temperature 
boundary layer profiles were found to be well-described 
by the Prandtl-Blasius profiles. Hence, our theory may 
explain the recent experimental observation of a reduc- 
tion in heat transport in turbulent RB convection with 
polymer additives. In particular, an amount of about 10 
% HR can be obtained with suitable parameters in our 
theory, and this amount of heat reduction is comparable 
to that observed in the experiment. 
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TABLE I: Amount of drag enhancement and heat reduction 
for different values of parameters. 



[1] K.R. Sreenivasan and C. White, J. Fluid Mech. 409, 149 
(2000). 

[2] I. Procaccia, V.S. L'vov, and R. Benzi, Rev. Mod. Phys. 

80, 225 (2008). 
[3] G. Ahlers and A. Nikolaenoko, Phys. Rev. Lett. 104, 

034503 (2010). 
[4] E.D. Siggia, Ann. Rev. Fluid Mech. 26, 137 (1994). 
[5] S. Grossmann and D. Lohse, J. Fluid Mech. 407, 27 

(2000). 

[6] Q. Zhou, R.J. A.M. Stevens, K. Sugiyama, S. Gross- 
mann, D. Lohse, and K.-Q. Xia, "Prandtl-Blasius tem- 
perature and velocity boundary layer profiles in turbulent 
Rayleigh-Benard Convection", J. Fluid Mech. 664, 297 
(2010). 

[7] See, for example, L.D. Landau and E.M. Lifshitz, Fluid 



Mechanics (Pergamon Press, Oxford, 1987). 

[8] H. Schlichting and K. Gersten, Boundary-Layer Theory 
(Springer, 8th ed. 2004). 

[9] R.B. Bird, O. Hassager, R.C. Armstrong, and C.F. Cur- 
tis, Dynamics of Polymeric Liquids (Wiley-Interscience, 
1987). 

[10] A. Celenti, A. Puliafito, and D. Vincenzi, Phys. Rev. 
Lett. 97, 118301 (2006). 

[11] D.O. Olagunju, App. Math. Lett. 19, 432 (2006). 

[12] R.C. Bataller, R.C, Phys. Lett. A 372, 2431 (2008). 

[13] R. Benzi, E.S.C. Ching, and V.W.S. Chu, "Polymer Ef- 
fects on Heat Transport in Laminar Boundary Layer 
Flow" , larXiv:1104.4230l 



